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Abstract 

We consider a special relativistic effect, known as the Poynting-Robertson effect, on various 
types of trajectories of solar sails. Since this effect occurs at order f'^/c, where v'^ is the 
tangential speed relative to the sun, it can dominate over other special relativistic effects, 
which occur at order v'^ je?. While solar radiation can be used to propel the solar sail, the 
absorbed portion of it also gives rise to a drag force in the tangential direction. For escape 
trajectories, this diminishes the cruising velocity, which can have a cumulative effect on the 
heliocentric distance. For a solar sail directly facing the sun in a bound orbit, the Poynting- 
Robertson effect decreases its orbital speed, thereby causing it to slowly spiral towards the 
sun. We also consider this effect for non-Keplerian orbits in which the solar sail is tilted in 
the azimuthal direction. While in principle the drag force could be counter-balanced by an 
extremely small tilt of the solar sail in the polar direction, periodic adjustments are more 
feasible. 



1 Introduction 



Solar electromagnetic radiation can be transmitted, reflected, absorbed and emitted by re- 
radiation by a solar sail. It is well known that the reflected, absorbed and emitted portions 
of the radiation can be used to propel the solar sail, due to the force of the electromagnetic 
pressure. What is less known is that the absorbed and emitted portions of the radiation 
induce a drag force on the solar sail, thereby diminishing its tangential speed relative to the 
sun. While this force is relatively small, it can have a long-term cumulative effect on the 
trajectories of solar sails. In particular, the effect of this drag force can be rather dramatic 
on the heliocentric distance of solar sails in long-range escape trajectories, while solar sails 
in bound orbits will spiral inwards to the sun. 

The aforementioned drag force on solar sails is associated with the Poynting-Robertson 
effect, which was flrst investigated for small spherical particles by Poynting in 1904 [1]. 
Although it is now realized that this is a special relativistic effect associated with the flnite 
speed of light, Poynting's paper was actually published a year before Einstein's paper on 
special relativity. This effect was reconsidered in 1937 by Robertson as a manifestly special 
relativistic effect |2]. Then starting from a paper by Wyatt and Whipple [3J, the Poynting- 
Robertson effect was used for modeling the orbital evolution of dust particles orbiting the 
sun and, in particular, accounted for the drag force which causes dust particles to slowly 
spiral inward towards the sun. Various reflnements of the original analysis have been made, 
such as including the leading-order effect of curved spacetime around the sun for particles 
of various shape^j. 

While there has been much research devoted to the influence of the Poynting-Robertson 
effect on the motion of dust grains, the Poynting-Robertson effect was only recently consid- 
ered for solar sails in bound orbits [6]. In this paper, we consider the Poynting-Robertson 
effect for solar sails in escape trajectories and further elaborate on the cases of bound he- 
liocentric and non-Keplerian orbits. From the rest frame of the solar sail, this effect is 
associated with the nonradial component of the absorbed solar radiation, due to the relative 
motion between the sun and the solar sail. From the rest frame of the sun, this results from 
the solar sail absorbing solar radiation and then emitting it in the forward direction relative 
to its motion. We would like to mention that there is a second effect that is also due to 
the absorption of solar radiation, known as the Yarkovsky effect which is related to the 
emissivity of the solar sail material. For a solar sail that has tangential motion, both types 
of effects can occur simultaneously. 

This paper is organized as follows. In section 2, we present the orbital equations for 
a solar sail orbiting in the plane of the sun, which takes into account the drag from the 

^See [4] and [5] and references therein. 
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Poynting-Robertson effect. In section 3, we investigate the diminished hehocentric distance 
which results from this drag effect on solar sails in escape trajectories. In section 4, we 
consider how the Poynting-Robertson effect causes an initially circular orbit to slowly spiral 
in towards the sun. We consider both bound orbits that lie within the plane of the center of 
mass of the sun, as well as non-Keplerian orbits which lie outside of this plane. Conclusions 
follow in section 5. 



2 Orbital equations 

The total force exerted on the solar sail due to solar radiation is the result of reflection, 
absorption and emission by re-radiation. When the solar electromagnetic radiation inter- 
acts with the solar sail material, it undergoes diffuse as well as specular reflection. The 
acceleration due to the diffuse reflection is directed along the normal to the antisun surface 
area, while the acceleration due to the specular reflection is directed opposite to the reflected 
radiation. In addition, the acceleration produced by the absorption of the solar radiation is 
directed along the incident radiation. This is shown in Figure [TJ Moreover, a portion of the 
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Figure 1: Specular reflection leads to an acceleration directed in the opposite direction from the reflected 
radiation, while absorption leads to an acceleration directed in the same direction as the incident radiation. 

absorbed radiation will be re-radiated from the front and back sides of the solar sail, which 
leads to an acceleration normal to whichever surface has the smaller emissivity. 

Following [SlIS], the acceleration due to the forces produced by the solar radiation pressure 
can be expressed as 

K L 
a = — cos^/'i [air - (2a3 cost/'i + a2)n] , k = ^ " , (1) 

Zuca 

where r is the unit vector in the direction of the heliocentric vector r, n is the unit normal 
vector for the surface area facing the sun, ipi is the pitch angle of the solar sail relative to the 
heliocentric vector r, a is the mass per surface area of the solar sail, and Ls= 3.842 x 10^^ 
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W is the solar luminosity at 1 AU. The optical characteristics of the solar sail film ai, a2, 03 
are defined as 



ai = 1 — ps, 



as = Bf{l-s)p+{l-p) 



e/ + efe 



(2) 



as = ps, 



where p is the reflection coefficient, s is the specular reflection factor, e/ and th are the 
emission coefficients of the front and back sides of the solar sail, and Bf and Bi, are the non- 
Lambertian coefficients of the front and back sides, which describe the angular distribution 
of the emitted and the diffusely reflected photons. For the case in which the surface of a 
non-perfectly reflecting solar sail is directly facing the sun, its acceleration is given by 



Equations ([T]) and ([3]) are well known. Now we will bring into consideration an additional 
force due to absorption which has not been taken into account in these equations. This is 
because the above description is in the Newtonian framework, for which the solar radiation is 
purely radial. Thus, for the case in which the surface of the solar sail is directly facing the sun, 
the force due to the solar radiation pressure (SRP) is directed radially outwards. However, 
as a special relativistic effect, the solar radiation has a nonzero tangential component in the 
frame of reference of the solar sail. This is due to the relative tangential speed v'^ = rcf) 
between the solar sail and the sun, where r is the heliocentric distance and is the angular 
coordinate. In particular, in the rest frame of the solar sail, the solar radiation propagates 
at an angle a = sin~^(f'^/c) with respect to the radial direction. Therefore, the absorbed 
portion of this radiation leads to a force with a component opposite the direction of motion. 
This drag effect is generally known as the Poynting- Robertson effect, which occurs at order 
v'^/c, and dominates over other special relativistic effects which are at order v'^/c^. Note that 
we are neglecting the redshift in wavelength of the solar radiation due to the radial velocity 
of the solar sail, which is a higher-order effect. 

We will first consider a solar sail whose surface is directly facing the sun and whose motion 
is restricted to lie within the heliocentric plane. The portion of light reflected by the solar 
sail leads to a radially outwards acceleration, whereas the portion of light absorbed leads to 
an acceleration directed at an angle a with respect to the radial direction, as depicted in 
Figure El 

We will be using the reflectivity parameter 0.5 < rj < 1 where rj = 0.5 (1) corresponds to 
the total absorption (total reflection) of solar radiation by the surface of the solar sail. The 
fraction of light reflected is 2?] — 1 and the fraction of light absorbed is 2(1 — 77). Throughout 
this paper, we will take a conservative value of // = 0.85. The acceleration due to the 




^ [air - (2a3 + a2)n] , 



(3) 
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Figure 2: Reflected light yields an acceleration in the radial direction, whereas absorbed light leads to an 
acceleration at an angle a with respect to the radial direction, due to the Poynting- Robertson effect. 

portion of light that is reflected by the solar sail is directed in the radial direction and has 
a magnitude of 

{2r]-l)K 2 

are fleeted = ^ COS a. (4) 

We are not considering fluctuations in the solar radiation and are taking a time-averaged 
luminosity. One factor of cos a is due to the fact that the area of the solar sail is effectively 
reduced since it lies at an angle of a with respect to the direction of the solar radiation. 
The second factor of cos a comes from projecting the solar radiation force vector along the 
direction normal to the surface of the solar sail. 

There is also acceleration due to the portion of light that is absorbed, which is directed 
at an angle a in the backward direction and has a magnitude of 



(1 — ri)Kcosa , , 

(^absorbed • y'^ ) 

Projecting the two contributions to the acceleration given by (j4]) and ([5]) along the radial 
and angular directions and including the acceleration due to the gravitational fleld of the 
sun gives 

rjK cos^ a GM 



r 

a 



[1 — r])^ sin a COS a , 
72 • (6) 
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where the mass of the sun is about M = 1.99 x 10 kg. 

In polar coordinates, the components of acceleration can be written as 



= r — rcj)^, 



a* = r-^(9t(r20). (7) 

Equating the components of acceleration in (|6]) and ([7]) yields 

GM — r]K cos^ a ■ ^ 
f H — r<p =0, 

9.(r^0) + ^i^^^l^^^^^^^ = O. (8) 
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Recalling that a = sin '^{v'^/c), we can write this up to leading order in v'^/c ^ 1 as 



GM-r]K -2 
r H — r(p =0, 

airV+^ = 0. 9 

cr 



Since v'^ = r(j), we can integrate the second equation to give 

• voro - ?7)(0 - 0o) 



(10) 



where ro, 0o and vq are the radial and angular positions and speed of the solar sail at closest 
approach to the sun. 

3 Escape trajectories 

We will take the initial conditions to be at closest approach, where the solar sail is deployed. 
Perihelion distances as small as 0.02 AU- 0.1 AU may be feasible for solar sails in the near 
future. For example, a trajectory design for a solar probe was presented in that includes 
repeated pole-to-pole sun flybys at a perihelion of slightly less than 0.02 AU. 

Before the solar sail is deployed at the distance of closest approach r = tq, the gravita- 
tional attraction of the sun causes the speed of the spaceship to increase as it gets closer to 
the sun. The Helios deep space probes would have traveled at the record speed of about 70 
km/s at 0.3 AU. This enables us to extrapolate that the following sampling of speeds fo are 
feasible for the near future: Vq = 133 km/s at tq = 0.1 AU, Vq = 188 km/s at tq = 0.05 AU 
and Vq = 298 km/s at ro = 0.02 AU. We will use these as sets of initial conditions in order 
to demonstrate the Poynting- Robertson effect, though of course our orbital equations can 
be applied to any initial heliocentric distance and velocity. We take the conservative value 
of 0" = 0.001 kg/m^ for the mass per area of the solar sail, since this could be large enough 
to take into account the mass of the load that is being transported. 

Most of the acceleration for a solar sail takes place during the first day after it has 
been deployed. For our three sets of initial conditions, the decrease in speed Av due to the 
Poynting- Robertson effect for the first day of the voyage is shown in the Figure [31 

For a solar sail deployed at 0.02 AU, the cruising speed after a 30-year voyage would be 
about 314 km/s. Throughout almost all of the voyage, the speed is about Av =20 m/s less 
than what it would have been in the absence of this effect, which then has a cumulative 
effect on the heliocentric distance. The Poynting-Robertson effect decreases the heliocentric 
distance by an amount of Ar 20 million kilometers after a 30-year voyage. For a solar 
sail deployed at 0.05 AU and 0.1 AU, the heliocentric distance is lessened by about 7 million 
kilometers and 4 million kilometers, respectively, compared to what it would have been in 
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Av (m/s) 




Figure 3: The decrease in speed Av due to Poynting- Robertson effect for a solar sail initially moving at 
0.02 AU with a tangential speed of 100 km/s (dashed line), 300 km/s (solid line) and 600 km/s (bold line). 

the absence of the Poynting- Robertson effect. Ar is shown as a function of the voyage time 
in Figure H] for our three sets of initial conditions. 

Ar (km) 




t (years) 



5 10 15 20 25 30 

Figure 4: The decrease in heliocentric distance Ar due to Poynting- Robertson effect versus the duration of 
the voyage for the following sets of initial conditions: vq — 133 km/s at tq = 0.1 AU (dashed line), vq = 188 
km/s at To = 0.05 AU (solid hne), and va = 298 km/s at tq = 0.02 AU (bold hne). 

4 Bound orbits in the plane of the sun 

We will now consider the Poynting-Robertson effect on a solar sail-propelled satellite in a 
bound orbit. As an example, we take the effective areal mass to be a = 0.00111 kg/m^. 
Table 1 lists the percentage decrease in the heliocentric distance after one year for a solar 
sail directly facing the sun and in a bound orbit at various initial distances from the sun. 
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Initial Speed 


Initial Distance 


Distance After 1 Year 


% Decrease 


781.17 m/s 


0.38 AU 


0.37986 AU 


0.04% 


1076.76 m/s 


0.2 AU 


0.1985 AU 


0.8% 


1522.77 m/s 


0.1 AU 


0.096 AU 


4% 


2153.52 m/s 


0.05 AU 


0.043 AU 


14% 


2407.71 m/s 


0.04 AU 


0.027 AU 


33% 


2780.19 m/s 


0.03 AU 


0.0044 AU 


85% 



The initial tangential velocity is taken to be the value necessary for having a circular orbit 
in the absence of the Poynting- Robertson effect. As can be seen, the amount of change in 
the heliocentric distance increases dramatically as the initial heliocentric distance decreases. 
For a solar sail-propelled satellite initially around the orbit of Mercury at about 0.38 AU, 
the heliocentric distance decreases by only 0.04% after one year. On the other hand, for an 
initial distance of 0.03 AU, the heliocentric distance decreases by about 85% after one year. 
In the extreme scenario of a satellite beginning at 0.02 AU, the heliocentric distance will 
decrease to less than 0.01 AU in about 0.7 years, as shown in Figure |5l 




-0.02 



Figure 5: Decrease in heliocentric distance r due to the Poynting-Robertson effect for a solar sail with 
(J = 0.00111 kg/m^ that is initially in a circular orbit at 0.02 AU. 

We will now consider a more general scenario in which the solar sail is tilted in the 
forward tangential direction (with respect to the velocity of the solar sail) at an angle ipi- 
Then the acceleration due to the portion of light that is reflected by the solar sail is given 
by 

(277 — 1)k 2/ , \ /i . N 

ttreflected = ^ COS (« + t/'i) , (11) 
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and is directed at an angle ipi in the forward direction. The acceleration due to the portion 
of light that is absorbed is directed at an angle a in the backward direction and has a 
magnitude of 

_ (1 -r/)/tcos(a + ^i) 

(^absorbed ^2 ' \ ) 

Projecting these two contributions to the acceleration along the tangential direction gives 
, [(2?7 — 1) cos(a + ipi) sinipi — (1 — r/) sina]fi;cos(a + ipi) 



(13) 



In principle, one could use the contribution to the acceleration in the forward tangential 
direction due to the reflected radiation to counteract the drag due to the absorbed radiation. 
In order to maintain a circular orbit at constant speed, we set a"^ = 0, which gives the 
condition 

(2?7 — 1) cos(a; + ipi) sin-i/'i = (1 — r/) sin a. (14) 

Since a = sin~^(t>/c), this relates the tilting angle ipi with the speed v and reflectivity 
parameter t] of the solar sail. For f ^ c, we can approximate ipi as 

For example, in order to cancel the drag for the solar sail with an initial velocity of 3405.02 
m/s (and t] = 0.85), we find that ipi ~ 2.4 x 10~^ rad. Since such precision for the tilting 
angle is clearly not realistic, one could instead perform periodic corrections on the tilting 
angle when the solar sail deviates too far from a circular orbit. 



5 Bound orbits out of the plane of the sun 

We will now consider non-Keplerian orbits, for which the center-of-mass of the sun does not 
lie within the orbital plane [IT1[T2]. We orient our coordinate system such that a circular 
non-Keplerian orbit would correspond to constant r and 6, with the orbital radius given by 
p = rsin6'. In the absence of the Poynting-Robertson effect, a non-Keplerian orbit would 
be maintained with a suitable pitch angle ip2 in the 6 direction and relative to the radial 
direction at the location of the solar sail, as shown in Figure |6l 

As before, due to the Poynting-Robertson effect, the solar radiation is at an angle a 
relative to the radial direction. We will consider a solar sail whose normal direction n to 
its surface is tilted at an angle ipi in the (j) direction, and which has a pitch angle ?/'2- We 
define \1/ as the angle between n and the radial direction, $ as the angle between n and the 
(f) direction, and 7 as the angle between n and the direction of the solar radiation. These 
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Figure 6: In the absence of the Poynting- Robertson effect, orienting a solar sail at an appropriate pitch 
angle i^j in the polar direction would lead to a circular orbit outside of the plane of the sun. 

angles are related to ^/^i, '^2 and a as follows: 



cos \1/ 
cos $ 
cos 7 



1 + tan^ t/^i + tan^ -^2 



a/1 + cot^ -^x + tan^ i\}2 ' 
1 

a/1 + tan^(a + '\\)]) + tan^ i\)2 



(16) 



The fraction of light reflected 21] — 1 pushes the solar sail at an angle relative to the 
radial direction, and the fraction of light absorbed 2(1 — rj) pushes the solar sail radially 
outwards. Thus, after performing the appropriate projections, we find 



GM 



[21] — 1)k cos^ 7 sin'?/'2 



[21] — l)/t cos^ 7 cos $ — (1 — cos 7 sin a 



(17) 



where 

M = M — — [(1 — 77) cos a + {21] — 1) cos 7 cos \E'] cos 7 . 
Acceleration in terms of spherical coordinates can generally be written as 



a 

a' 

at' 



r — rO'^ — r sin^ 9 0^ , 
r9 + 2r6 — r sin cos di^ 



sin 9 



dt(r'^(j)) + 2r cos 6' 



(19) 



where ' = dt. Equating the components of acceleration in (fTTj) and f|T9l) yields the equations 
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of motion for the solar sail: 

r ~ rU — r sm 6 (p = — , 



2 , 



a , o-a -a a (2r/ - l)Kcos2 7sin?/^2 
ro + 2r6 — r sm 6 cos a (p = , (20) 



r 



sin6'_,n-, „ ^ •• (2r7 — l)/tcos^7COS$ — (1 — ?7)fi:cos7sina 
dtir^d)) + 2r cos e ^-'^ _ v ' ^ v /y 
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Figure 7: A three-dimensional trajectory for a solar sail initially in a circular orbit outside of the plane of 
the sun at 0.02 AU and a polar angle of 88.28°. 

Figure [7] shows a three-dimensional trajectory during one year, due to the Poynting- 
Robertson effect on a solar sail with a = 0.00111 kg/m^ that is initially in a non-Keplerian 
orbit at tq = 0.02 AU and = 88.28° with an initial speed of Vq = 3406.63 m/s and a pitch 
angle of ip = 0.0054°. In order to depict the different types of motion that can occur in such 
a trajectory, Figure [8] shows an example trajectory along the r, 6 and directions for which 
the solar sail is initially at tq = 0.02 AU and 6q = 84.27° with an initial speed of Vq = 3422.12 
m/s and a pitch angle of ip = 0.0018°. As can be seen, the heliocentric distance r is reduced 




Figure 8: Spherical coordinates of a three-dimensional trajectory for a solar sail initially in a circular orbit 
outside of the plane of the sun at 0.02 AU and a polar angle of 84.27°. 

to 0.01 AU in 0.7 years. Since the solar sail gets closer to the sun, the azimuthal velocity in 
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the (f) direction increases. Notice also that the solar sail exhibits oscillatory behavior in the 
9 direction, which can be minimized by tuning the pitch angle. 

6 Conclusions 

We have considered the Poynting- Robertson effect resulting from the absorption of the solar 
radiation by a solar sail. The Poynting- Robertson effect occurs at order v'^/ c, and dominates 
over other special relativistic effects which are at order v'^/c^. For escape trajectories, this 
effect decreases the cruising velocity as well as the heliocentric distance. For the example 
of a solar sail with mass per area a — 0.001 kg/m^ deployed between 0.02-0.1 AU, the 
Poynting-Robertson effect lessens the hehocentric distance after a 30-year voyage by 4-20 
million kilometers. 

As for bound orbits, if not compensated for, this effect can decrease the orbital speed of 
the solar sail, thereby causing it to slowly spiral towards the sun. We considered an example 
of a solar sail whose surface was directly facing the sun with mass per area a = 0.00111 
kg/m^. The percentage by which the heliocentric distance changes due to the Poynting- 
Robertson effect increases dramatically as the initial distance decreases. For an initial dis- 
tance of 0.1 AU, the heliocentric distance decreases by about 4% in one year, whereas an 
initial distance of 0.03 AU is decreased by 85% in one year. In principle, the drag effect 
could be counter-balanced by an extremely small tilt of the solar sail in the polar direction. 
However, periodic adjustments when the solar sail deviates too much from the circular orbit 
are more feasible. 

We have also considered the Poynting-Robertson effect for a solar sail with a three- 
dimensional orbit. As an example, we have considered a solar sail with the necessary values 
of parameters so that, in the absence of the Poynting-Robertson effect, it would follow a 
closed circular orbit above the plane of the sun at a hehocentric distance of 0.02 AU. While 
it spirals closer towards the sun due to the Poynting-Robertson effect, it also undergoes 
oscillatory motion in the polar direction, which can be controlled by the pitch angle. 

For simplicity, we have worked in the approximation that the sun is a point-like light 
source. One could also take into account the angular variations of the solar radiation asso- 
ciated with the shape of the sun. However, this will not change our conclusion about the 
importance of the Poynting-Robertson effect. 

Table 1: The percentage decrease in hehocentric distance due to the Poynting-Robertson 
effect for bound orbits with various initial conditions. 
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